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ABSTRACT 

We present Spitzer 7-38 /im spectra for a 24 pm flux limited sample of galaxies at z ~ 0.7 in the COSMOS 
field. The detailed high-quality spectra allow us to cleanly separate star formation (SF) and active galactic 
nucleus (AGN) in individual galaxies. We first decompose mid-infrared Luminosity Functions (LFs). We find 
that the SF 8 /im and 15 /im LFs are well described by Schechter functions. AGNs dominate the space den- 
sity at high luminosities, which leads to the shallow bright-end slope of the overall mid-infrared LFs. The 
total infrared (8-1000 /im) LF from 70 /im selected galaxies shows a shallower bright-end slope than the bolo- 
metrically corrected SF 15 /im LF, owing to the intrinsic dispersion in the mid-to-far-infrared spectral energy 
distributions. We then study the contemporary growth of galaxies and their supermassive black holes (BHs). 
Seven of the 3 1 Luminous Infrared Galaxies with Spitzer spectra host luminous AGNs, implying an AGN duty 
cycle of 23 ± 9%. The time-averaged ratio of BH accretion rate and SF rate matches the local MBH-A^buige 
relation and the MBH-A^host relation at z ~ 1 . These results favor co-evolution scenarios in which BH growth 
and intense SF happen in the same event but the former spans a shorter lifetime than the latter. Finally, we 
compare our mid-infrared spectroscopic selection with other AGN identification methods and discuss candi- 
date Compton-thick AGNs in the sample. While only half of the mid-infrared spectroscopically selected AGNs 
are detected in X-ray, —90% of them can be identified with their near-infrared spectral indices. 
Subject headings: galaxies: active — galaxies: luminosity function — quasars: general — infrared: galaxies 
— X-rays: galaxies 



1. INTRODUCTION 

A central question in extragalactic research is the mass as- 
sembly history of galaxies and their supermassive black holes 
(BHs). Understanding this process requires measurements of 
the rates of star formation (SF) and BH accretion throughout 
the cosmic ages. As both processes typically occur in dusty 
environments, the majority of the energy emerges as dust- 
reprocessed thermal infrared (IR) emission. The IR luminos- 
ity function (LF), defined as the density distribution over IR 
luminosity, thus holds clues to the evolution of galaxies and 
BHs. 

Measuring IR LFs has a long history. About two decades 
ago, surveys with the Infrared Astronomical Satellite (IRAS) 
established the local benchmarks in the mid-IR and far-IR 
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wavelengths (A > 5 pm; e.g., Rush et al. 1993; Soifer et al. 
1987). Although the space density of Luminous and Ultra- 
Luminous Infrared Galaxies (LIRGs/ULIRGs; log(LiR/L Q ) = 
log[L(8-1000/im)/L Q ] = [11, 12]/[12, oo)) exceeds those of 
optically-selected galaxies and QSOs at comparable bolomet- 
ric luminosities, they account for only —5% of the total in- 
tegrated IR energy density in the local universe (Soifer & 
Neugebauer 1991). LIRGs, however, dominate SF activities 
at higher redshifts. Strong evolution of IR-selected popula- 
tion with look-back time was suggested by the number count 
results from the Infrared Space Observatory (ISO) observa- 
tions (Elbaz et al. 1999). Remarkably, this evolution was first 
detected in the limited redshift range covered by the earlier 
IRAS surveys (0 < z < 0.2; e.g., Saunders et al. 1990). Re- 
cently, many attempts have been made to estimate IR LFs at 
redshifts up to z ~ 3 using data from the new generation of 
IR space telescopes — Spitzer (Le Floc'h et al. 2005; Perez- 
Gonzalez et al. 2005; Franceschini et al. 2006; Bell et al. 
2007; Caputi et al. 2007; Magnelli et al. 2009; Rodighiero 
et al. 2010), AKARI (Goto et al. 2010), and Herschel (Grup- 
pioni et al. 2010; Eales et al. 2010). The strong evolution was 
confirmed; Le Floc'h et al. (2005) concluded that the 24 fim- 
derived IR (8-1000 /im) comoving energy density evolves as 
(1 +z) 4 at < z < 1 and LIRGs contribute —70% of the energy 
density at z = 1 . 

Although a crude picture has already emerged, our knowl- 
edge of the galaxy LFs is still limited by two main fac- 
tors: (1) the poorly sampled IR spectral energy distribution 
(SED), and (2) the elusive contribution from active galactic 
nuclei (AGNs). Constrained by the low sensitivities of ex- 
isting far-IR and sub-milimeter instruments, previously pub- 
lished 8-1000 )ira LFs at z > 0.3 are mostly /^-corrected mid- 
IR LFs. The /f-correction depends strongly on the assumed 
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SED, which without long wavelength data one can only guess 
from the SEDs of local galaxies. By opening the far-IR win- 
dow, Herschel will soon revolutionize this field. On the other 
hand, besides dust-enshrouded SF, mid-IR emission can also 
be powered by AGNs, as their accretion disks can heat the 
surrounding dusty tori. One therefore should separate AGN 
and SF before using mid-IR data to build 8-1000 pm LFs. 
X-ray imaging has been widely used to identify AGNs but 
the identification becomes incomplete at moderately high red- 
shifts, because of the insufficient depth in hard X-ray, which 
has less biases against absorbed AGNs. In addition, the sig- 
nificant population of AGN/SF composite systems renders it 
problematic to simply remove all of the identified AGNs from 
the sample. Only with mid-IR spectroscopy we can decom- 
pose AGN and SF in the mid-IR LFs by separating the two in 
individual sources. 

AGNs exhibit a class of IR SEDs distinct from star-forming 
galaxies (Hao et al. 2007; Richards et al. 2006; Netzer et al. 
2007; Fu & Stockton 2009), due to the clumpy nature of the 
torus (Sirocky et al. 2008; Nikutta et al. 2009). Broadly speak- 
ing, AGN SEDs (vLv) are flat from near-IR to mid-IR but 
drop beyond 20 /im, and SF SEDs peak between 55-120 /im 
(Rieke et al. 2009). 

The shape of the LF depends on the rest-frame wavelength. 
From UV to near-IR wavelengths, the classical Schechter 
(1976) function provides excellent fits (e.g., Arnouts et al. 
2005; Ilbert et al. 2005; Babbedge et al. 2006; Dai et al. 2009): 
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At mid-IR and far-IR wavelengths, the observed bright-end 
slope {<f> oc L~ 2 ) is shallower than the exponential cut-off of 
the Schechter function (e.g., Rush et al. 1993; Lawrence et al. 
1986). These LFs are best described by modified Schechter 
functions (Saunders et al. 1990): 
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It may appear counter-intuitive that UV LFs show a profile 
different from that of far-IR LFs, since both UV and far-IR 
emission trace active SF. One explanation is that since the UV 
LF is not extinction corrected, its bright end is steepened be- 
cause dust extinction increases with star formation rate (SFR; 
Hopkins et al. 2001; Martin et al. 2005; Zheng et al. 2007b; 
Reddy et al. 2010). 

AGNs might be responsible for the shallow bright-end slope 
of the mid-IR LFs. Rush et al. (1993) showed that the local 
12 /im LF of normal galaxies (i.e., excluding Seyferts of both 
spectral types) has a steeper bright-end slope (0 oc L~ 3 6 ). Ad- 
ditionally, after removing spectroscopically classified AGNs, 
Huang et al. (2007) found the local 8 /im LF was well fit by 
a Schechter function (see also Dai et al. 2009). At higher 
redshifts, the results are less conclusive, because of the limi- 
tations from both the small survey areas and the uncertainties 
in AGN identification. 

In this paper, we decompose AGN and SF in mid-IR LFs 
at z ~ 0.7 using a unique spectroscopic data-set of LIRGs. 
We extend these results to lower luminosities (Lir > 10 10 L©) 
with the 2 deg 2 S-COSMOS 24 pm catalog (Sanders et al. 
2007; Le Floc'h et al. 2009). In § 2 we present Spitzer In- 



frared Spectrograph (IRS; Houck et al. 2004) spectra of a 
24 /tm flux limited sample at z ~ 0.7 and cross-match COS- 
MOS multi-wavelength photometry and redshift catalogs. In 
§ 3, we model the IRS spectra to divide the observed lumi- 
nosities between AGN and SF. This method, for the first time, 
allows us to consider all of the sources in the sample and de- 
rive AGN and SF LFs separately (§ 4). To estimate the total 
SFR function, we use a 70 /im selected sample to derive the 
8-1000 /im LF and compare it with the SF mid-IR LF (§ 4.3). 
In § 5, we explore galaxy-BH co-evolution in LIRGs, com- 
pare our mid-IR spectroscopic AGN identification with pho- 
tometric AGN selections, and discuss the possible presence of 
Compton-thick AGNs in our sample. We close by summariz- 
ing our conclusions in § 6. 

Throughout we assume a cosmological model with Hq = 70 
km s" 1 Mpc" 1 , Q, m = 0.3, and £1a = 0.7. All magnitudes 
are in the AB system. The SFRs are given for the Chabrier 
(2003) initial mass function (IMF); they can be converted to 
the Salpeter (1955) IMF by adding —0.24 dex to the loga- 
rithm. 

2. DATA 

The data required to build mid-IR LFs include: (1) IRS 
spectroscopy of a complete 24 /im flux limited sample {F2A > 
0.7 mJy), (2) the Multiband Imaging Photometer for Spitzer 
(MIPS; Rieke et al. 2004) 24 /im catalog down to 60 /Jy 
(Sanders et al. 2007; Le Floc'h et al. 2009), and (3) the 
redshift catalogs (both spectroscopic and photometric red- 
shifts). The IRS spectra were used to decompose AGN and 
SF. The MIPS catalog was used to build LFs down to the 
luminosity limit of the S-COSMOS survey. In order to re- 
move stellar photospheric emission in the mid-IR, we also in- 
cluded the 30-band COSMOS photometry from Galex/Ft/V 
(155 1A) to Spitzer/IRAC 8 /im. For the 8-1000 /im LF, we 
used the MIPS 70/160 /im data (Frayer et al. 2009). In this 
section, we describe our IRS observations and explain the 
multi-wavelength cross-matching method. 

2. 1 . IRS Spectroscopy 

We selected the targets of our Cycle 5 SpitzerflRS Legacy 
program (PID 50286, PI: Scoville) using the following crite- 
ria: (1) they are inside the COSMOS field observed by Hubble 
Space Telescope/ 'Advanced Camera for Surveys (HST/ACS), 
(2) have MIPS 24 pm flux brighter than 0.7 mJy, and (3) 
have photometric redshifts between 0.6 < z < 0.8. Nine 
blended sources were removed to avoid confusion in the spec- 
tra. AORs were generated for 55 galaxies. 53 were success- 
fully observed, and 2 targets failed due to errors in the peak-up 
imaging acquisition (Table 1). 

At the time of writing, 48 of the 53 galaxies (90%) have 
spectroscopic redshifts (including eight new redshifts from 
the IRS spectra). In Figure 1 we show the distributions of 
this sample in the COSMOS field, in redshift, and in 24 pm 
flux. We compare these distributions with the parent sample 
using the currently best-estimated redshifts (§ 2.2). We have 
spectroscopic redshifts for 49% of the 24 pm sources above 
0.7 mJy and for 79% of the X-ray detected 24 /im sources 
above 0.7 mJy. Most of the spectra come from zCOSMOS 
(Lilly et al. 2007) and the COSMOS AGN Spectroscopic Sur- 
vey (Trump et al. 2009). We found 70 sources above 0.7 mJy 
between 0.6 < z < 0.8 over the 1.63 deg 2 HST/ACS field 
(Scoville et al. 2007). 48 of the 70 (68.6%) have IRS spec- 
tra. Accounting for this difference, the effective area of the 
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FIG. 1. — IRS sample selection and completeness, (a) Distribution of the IRS sample in the COSMOS field. The dashed polygon delineates the HST/ ACS 
coverage. Black open circles mark the 70 MIPS sources at 0.6 < z < 0.8 and F24 > 0.7 mjy inside the ACS field. 48 of those have IRS spectra and are filled in 
blue. Red open squares mark the 5 IRS sources that fall outside of the redshift window. The red crosses show the 2 IRS sources for which the observations failed 
due to peak-up acquisition errors, (b) Redshift distribution of the MIPS 24 fim sample with F24 > 0.7 mJy and inside the ACS field (black solid histogram), in 
comparison with that of the IRS sample (blue/red filled histogram). Grey shaded region shows our redshift selection window, (c) 24 fim flux vs. redshift. MIPS 
sources inside the ACS field with F24 > 0.7 mJy and F24 < 0.7 mJy are shown as black and grey points, respectively. The blue open circles and red open squares 
highlight the 53 sources with IRS spectra. 



IRS sample is 1.12 deg 2 . We included the 5 IRS sources out- 
side of 0.6 < z < 0.8 in the spectral analysis (§ 3), although 
we ignored them when estimating LFs. 

Low-resolution IRS spectra (R ^ 60- 130) were taken in 
January 2009. We used the Short Low module in 1st order 
(SL1, 7.4-14.5 fim) and the Long Low module in 1st (LL1, 
19.5-38 /im) and 2nd order (LL2, 14.0-21.3 fim) to contin- 
uously cover a wavelength range of 7.4-38 /im. All obser- 
vations were made in the standard staring mode. One cycle 
consisted of two nod positions offset by 20" and 55" for SL 
and LL modules, respectively. The number of cycles was 
adjusted to yield approximately equal signal-to-noise ratios 
across modules and sources. Between 2 and 25 cycles of 60 s, 
2 and 15 cycles of 120 s, 2 and 10 cycles of 120 s were used 
in SL1, LL2, and LL1, respectively (Table 1). 

The data were first passed through the S18.7 version of the 
SSC pipeline. Spectra were then extracted from the BCD 
files using the IRSLow IDL pipeline (Fadda et al. 2010). 
Briefly, the pipeline procedure is as follows: (1) After re- 
jecting bad pixels and subtracting background from individual 
frames, the pipeline combines all of the 2D spectra at differ- 
ent nod positions; (2) Optimal background subtraction is done 
in the second pass by manually masking off objects from the 
background-subtracted coadded image; (3) ID spectra are ex- 
tracted using PSF (Point Spread Function) fitting with PSF 
profiles from IRS calibration stars. Lastly, we removed bad 
portions near the edges of each order from the ID spectra. 

The LL spectra of MIPS 3 1635 and 46960 appear unusable. 
The former is contaminated by a nearby unidentified bright 
object (likely an asteroid), and the latter is dominated by a 
featureless continuum ~ 5 x higher than the MIPS 24 /im flux. 
We discarded these spectra and used their MIPS 24 /im pho- 
tometry as surrogates. 

In order to check the flux calibration, we measured 24 fim 
fluxes from the IRS spectra by convolving them with the 
MIPS filter transmission curve. Note that at z ~ 0.7, our ob- 
jects are unresolved in the 24 /im image, and the flux loss 



due to the finite slit widths have been fully accounted for 
during the spectral extraction. We found F24irs /F24mips = 
1 .078 ± 0.096. The MIPS fluxes are systematically lower than 
the IRS fluxes mainly because they are calibrated against a 
10,000 K blackbody instead of a flat F v as in the AB system. 
As we determined our ^-corrections for the 24 fim sources 
using the IRS spectra (§ 4.1), we multiplied the MIPS 24 fim 
fluxes by 1.078 for sources between 0.6 < z < 0.8. 

2.2. Multi-Wavelength Dataset 

Redshift is critical for estimating LFs. We matched the 
MIPS 24 fim catalog (Le Floc'h et al. 2009) with the COS- 
MOS dataset to identify redshifts. The full MIPS cata- 
log contains 52,092 sources down to 60 fiiy. The cross- 
matching was performed inside the 2 deg 2 Subaru deep area 
(149.4114075° < a < 150.8269348° and 1.4987870° < 5 < 
2.9127350°). The optical counterparts of 24 fim sources were 
identified in two steps. First, the 24 fim coordinates were 
matched to the nearest IRAC detection within a 2" radius 
(Ilbert et al. 2010). Then, we searched for the nearest opti- 
cal counterpart at < 26.5 within 1" to the IRAC position. 
We used the combined FUV-Xo-Ks photometric catalog (Ca- 
pak et al. 2007) so that the FUV-to-24 fim SEDs were assem- 
bled at the same time. When spectroscopic redshifts were un- 
available 13 , we used the photometric redshifts of Ilbert et al. 
(2009) and Salvato et al. (2009) for normal galaxies and X- 
ray sources, respectively. Lastly, we removed sources inside 
the Subaru/optical and SpitzerflRAC 3.6 fim masks, where the 
quality of photometry and photometric redshifts degrades due 
to the proximity to very bright objects. The area of this final 
catalog is 1.66 deg 2 (the masked area is 0.28 deg 2 ), enclosing 
a total of 39,525 MIPS sources down to 60 /iJy. Among these, 
34,967 have secure optical counterparts. 

12 Slit widths are 3."7, 10"5, and 10'.' 7 for SL1, LL2, and LL1 modules, 
respectively. 

13 15% of the 24 /im sources in the 60 /.tjy catalog have spectroscopic red- 
shifts. 
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FIG. 2. — The optical identification of 24 /im sources above 60 /Jy and between 0.6 < z < 0.8 is ~99% complete, (a) Identification completeness of 24 fim 
sources at all redshifts as a function of 24 fim AB magnitudes (solid line). Dotted histogram shows the distribution of the 39,441 objects within the Subaru 
deep area but outside of the masked areas, and the light grey filled histogram shows the distribution of the 34,883 objects with secure optical counterparts. The 
identification completeness is simply the ratio of the two. (b) Redshift distributions of 24 fim sources with optical counterparts at 20.5 < < 22.5 (black), 
22.5 < i'^ B < 24.5 (grey), 24.5 < i^ < 26.5 (light grey), and i^ < 26.5 (dotted histogram). 



We found 84 instances of source fragmentation in the MIPS 
catalog, and we combined the MIPS fluxes for these sources. 
The MlPS-selected sample was thus reduced to a total of 
39,441 objects with 34,883 secure counterparts (88.4%). Fig- 
ure 2a shows the optical identification completeness as a func- 
tion of 24 ^m magnitudes. Note that the completeness reaches 
75% even for objects as faint as [24]ab = 19.45 (= 60 /iJy). 

To build the integrated 8-1 000 fim LF, we also identified 
closest MIPS 70/160 /im counterparts of 24 fim sources us- 
ing the MIPS -Germanium 3a (peak-to-noise) catalog (Frayer 
et al. 2009). We adopted a matching radius that equals the 
cataloged 2a radial positional error from PSF fitting. Fol- 
lowing Kartaltepe et al. (2010), 559 spurious 70 fim sources 
were removed from the catalog. We found 70/160 fim coun- 
terparts for 1831/528 of the 34,883 24 fim sources that are 
optically identified. Our counterpart identification covers a 
^50% larger area than that of Kartaltepe et al. (2010), since 
we did not restrict it to the ACS field; but both 70 fim samples 
produce consistent 8-1000 ^m LFs (§ 4.3). 

3. DECOMPOSITION OF IRS SPECTRA 

Mid-IR spectra can be decomposed into SF and AGN com- 
ponents (e.g., Sajina et al. 2007; Pope et al. 2008). Photo- 
spheric emission from low-mass stars can be significant below 
10 fim. To estimate its contribution in the IRS spectra, we fit 
the rest-frame 912 A to 1.3 fim SED (i.e., GALEX FUV to K s 
bands) using Bruzual & Chariot (2003) stellar population syn- 
thesis models and extrapolated the best-fit SED to the mid-IR. 
Once the photospheric emission was subtracted, we modeled 
the IRS spectrum as a superposition of dust-obscured SF and 
AGN. The former is dominated by polycyclic aromatic hydro- 
carbons (PAH) features, and the latter can be approximated 
by a power law. For the SF component, we used the aver- 
aged templates of Rieke et al. (2009), built from nearby star- 
forming galaxies in 14 luminosity bins, 9.75 < log(LiR/L Q ) < 
13. Both the AGN and SF components were subjected to at- 
tenuation from a screen of dust, adopting the modified Galac- 
tic center extinction curve of Smith et al. (2007) with j3 = 0.1. 

Our fitting procedure is as follows. First, we corrected 
the templates for silicate extinction and stellar emission with 



PAHFit (Smith et al. 2007). We used the same extinction 
curve and approximated the stellar emission with a 5,000 K 
blackbody. Then, we found the best solution for each of the 

14 templates using MPFIT, an IDL ^-minimization routine 
(Markwardt 2009). As we allowed different amount of extinc- 
tions for the two components, there are five free parameters in 
each model — power-law index, scaling factors for the two 
components, and their extinctions at 9.7 fim (r^ N , t|^). Fi- 
nally, the solution giving the global minimum \ 2 was selected 
as the best two-component model. 

Since many sources appear dominated by either SF or AGN, 
we also performed least-% 2 fits with a single component. In 
the AGN-only model we fit for both silicate emission and ab- 
sorption. The x~ values from the two single-component mod- 
els were compared with that of the two-component model, 
and the best model was determined using an F test. The two- 
component model was deemed necessary only when its % 2 
reduction relative to the single-component models yielded an 
F value that is less than 5% likely (i.e., the probability that the 
X 2 reduction is due to random error is less than 5%). We used 
SF-only, AGN-only, and composite models in 22, 10, and 21 
sources, respectively. 

We sorted our modeling results by ascending rest-frame 

15 fim luminosity (Z45) for the 48 sources between 0.6 < z < 
0.8 in Figure 3. The AGN contribution increases dramati- 
cally with luminosity: while they are absent in less luminous 
sources (log(Li 5 ) < 10.6 L or F24 < 0.8 mJy), AGNs dom- 
inate above log(Li5) > 10.9 L Q (F24 > 1 .2 mJy). The sources 
hosting AGNs show IRAC fluxes clearly above the best-fit 
stellar population synthesis models; and the amount of the ex- 
cess agrees well with the best-fit power-law model. At longer 
wavelengths, since a typical AGN SED declines sharply be- 
yond rest-frame 20 fim, one would expect the far-IR emission 
to be dominated by SF. For the 15 MIPS 70 /xm-detected SF- 
only sources, although the extrapolation from the best-fit SF 
template generally over-predicts the 70 fim flux, the differ- 
ence is within 0.5 dex. The disagreement is more severe for 
the SF/AGN composite systems, and the three 70 /im-detected 
AGN-only sources all have 70 fim fluxes in excess of those ex- 
trapolated from the mean AGN SED of Netzer et al. (2007), 
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FIG. 3. — Decomposition of AGN and SF with IRS spectra. Panels are sorted by ascending rest-frame 15 fim luminosity to show the dramatic increase of AGN 
contribution with luminosity. Labeled are log(Li5/Lo) and MIPS IDs. IRS spectra are shown in grey and photometric data are green horizontal bars. The black 
curve shows our best-fit model to the IRS spectra, which is a combination of stellar photospheric emission (purple), dust-reprocessed SF (blue) and power-law 
AGN (red). 
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FIG. 4. — PAHFit 11.3 iim PAH EQW vs. rest-frame 15/im luminosity 
(Ll 5 ). The IRS-selected AGNs (color-filled symbols) have lower PAH EQWs 
and higher luminosity, similar to local warm ULIRGs (Genzel et al. 1998; 
Armus et al. 2007). X-ray detected sources are plotted as stars. The colors 
of the symbols are coded according to the AGN contribution integrated over 
5.5 < A re st < 20 fim (pagn), and the sizes increase logarithmically with the 
AGN luminosity at rest-frame 15 fim (L*™). The median errors are shown 
at the top left corner. The upper limits indicate spurious PAH detections. The 
solid and dotted lines show the mean and 1<t dispersion of the EQWs of pure 
star-forming galaxies (/?agn = 0). 



hinting that AGNs could also be a significant power source of 
the rest-frame 40 /mi emission. We will discuss further on the 
far-IR emission in § 4.3. 

As a measure of the strength of SF relative to that of AGN 
heated dusts, PAH equivalent width (EQW) is recognized as a 
good AGN diagnostic (Genzel et al. 1998; Armus et al. 2007). 
In Fig. 4 we show that the 11.3 /im PAH EQW decreases as 
AGN intensifies, confirming that our spectral decomposition 
and the PAH EQWs are equivalent in identifying AGNs. Note 
that the PAH EQWs were measured with PAHFit (Smith et al. 
2007), which yields EQWs a few times larger than the spline 
method (Fadda et al. 2010). We will present comparisons with 
other commonly used AGN diagnostics in § 5.2. 

4. LUMINOSITY FUNCTIONS 

4. 1 . Completeness, K-correction and lA/ msa Estimator 

Before estimating the LFs, one has to take into account 
three major sample biases: (1) source extraction incomplete- 
ness, (2) optical identification or redshift incompleteness, and 
(3) photometric redshift uncertainties. Monte Carlo simula- 
tions showed that source extraction from the 24 /im mosaic 
is highly complete; the catalog is about 75% complete at 
F2A = 60 /iJy and reaches 100% completeness above 1 10 /iJy 
(Le Floc'h et al. 2009). The completeness of the 70 /im cat- 
alog is poor; it is only ^20% complete at 7 mJy and ap- 
proaches 100% completeness above 16 mJy (Frayer et al. 
2009). The optically unidentified fraction of 24 /im sources at 
z ~ 0.7 should be less than 1%, since only 0.6% of the MIPS 
sources between 0.6 < z < 0.8 have optical counterparts be- 
tween 24 < /+ B < 26.5 (Fig. 2b; Le Floc'h et al. 2009). The 
photometric redshift catalogs show excellent agreement with 
spectroscopic redshifts for sources with z^ B < 24 and z < 1 .25 
— <ta.-/(i+;) = 0.012(0.015) for normal galaxies and X-ray 
sources, respectively (Ilbert et al. 2009; Salvato et al. 2009). 
Since 99.4% of the 24 /im sources at 0.6 < z < 0.8 have op- 
tical counterparts brighter than < 24, photometric redshift 



scattering should not affect our results. Besides, 29.4% of the 
24 /tm sources down to 60/iJy between 0.6 < z < 0.8 have 
secure spectroscopic redshifts. To conclude, the only sample 
bias that we need to correct for is from source extraction. 

With the IRS spectra, we can compute LFs at any wave- 
length between 5.5 and 20 /im. We decided to show LFs at 
two representative mid-IR wavelengths, 15 /im (monochro- 
matic) and IRAC 8 /im, for several reasons: (1) 24 /im roughly 
corresponds to rest-frame 15 /im at z ~ 0.7, minimizing band- 
pass correction, (2) for star-forming galaxies the 15 /im flux is 
dominated by dust continuum from very small grains stochas- 
tically heated by young stars and the IRAC 8 /tm filter in- 
cludes strong PAH features at 7.7 and 8.6 /im, and (3) local 
benchmarks at these wavelengths exist in the literature (e.g., 
Xu 2000; Pozzi et al. 2004; Huang et al. 2007). 

As the luminosity limit of the IRS sample is a few times 
higher than the predicted "knee" of the LF, it is necessary to 
include the entire 24 /im sample to complement these bright- 
end LFs. Without mid-IR spectra for these sources, we have 
to ^-correct the observed 24 /im fluxes to rest-frame lumi- 
nosities using assumed SEDs. Since AGNs are less relevant 
at lower luminosities (Figs. 3 & 4), it is reasonable to assume 
that sources below 0.7 mJy are dominated by star-forming 
galaxies (Fadda et al. 2010). Hence, we used the IRS spectra 
of SF-dominated galaxies to select the best SED for the K- 
correction. We compared the measured /^-corrections of the 
20 "pure" star-forming galaxies (i.e., mid-IR AGN contribu- 
tion, /?agn> l ess man 8%) with predictions of the templates in 
four commonly used SED libraries for IR galaxies (Chary & 
Elbaz 2001; Dale & Helou 2002; Lagache et al. 2003; Rieke 
et al. 2009). In Figure 5, we show, for each SED library, the 
track of the template that yielded the smallest residual to the 
average measured /T-correction at each redshift. Although we 
selected the best fit of the four best templates 14 for subse- 
quently ^-correcting the 24 /tm sources, we obtained almost 
identical results for all four of them. 

We adopted a single template for the ^-correction because 
of practical reasons: (1) no trend of ^-corrections relative to 
luminosity can be detected because of the narrow luminosity 
range spanned by the pure star-forming galaxies (logtLis/L©) 
= 10.6-11.1), and (2) templates that match the observed 
[8.0]-[24] colors at a given redshift produce erroneous K- 
corrections at 8 and 15 /tm. 

We computed LFs with the \IV mdx method (Schmidt 1968). 
We calculated the lower and upper redshift limits {znan and 
Zmax) for each galaxy to remain above the MIPS flux limits. 
Since MIPS sources at z <~ 0.7 are bright in the optical, the 
effect of the optical selection (/^ B < 26.5) is negligible. We 
applied weights (defined as 1/completeness) in the \IV max cal- 
culation to correct for sample incompleteness below 0.11 and 
16 mJy for 24 and 70 /tm, respectively. Errors of the comov- 
ing densities were calculated based on Poisson statistics. 

4.2. Mid-IR Luminosity Functions 

The total 15 /tm LF from the MIPS sample agrees with that 
of Le Floc'h et al. (2005). It is consistent with a modified 
Schechter profile (Eq. [2]), although our finer luminosity bins 
revealed a discontinuity at log(Li5/L Q ) ~ 10.8. Accounting 
for the different effective areas, the LFs from the MIPS sam- 
ple and the IRS sample agree well with each other, confirming 
that the IRS sample is representative of the complete sample. 

14 A Rieke et al. (2009) template (L iR = 5.3 X 10 10 L ) at 15 fim and a 
Lagache et al. (2003) template (Lir= 9.1 X 10 12 Lq) at 8 /im. 
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FIG. 5. — (a) 15 /im and (i) IRAC 8 fim /T-corrections from MIPS 24 fim. The black points show A"-corrections of the IRS sample as measured from their 
mid-IR spectra. The red boxes highlight the 20 pure star-forming galaxies (/>agn = 0). For each of the four IR libraries (Lagache et al. 2003; Rieke et al. 2009; 
Dale & Helou 2002; Chary & Elbaz 2001), we show the /f-correction track of the template that best fit the measured A'-corrections from the 20 star-forming 
galaxies. The selected templates for X'-corrections are highlighted with thicker lines. 



Using our mid-IR decomposition results (§ 3), we measured 
AGN and SF luminosities separately for individual sources. 
The SF and AGN LFs are shown in Fig. 6 and are tabulated in 
Tables 2 & 3. 

The AGN subtraction removed the relatively large number 
of high-luminosity sources, which had justified the usage of 
modified Schechter profiles for the overall mid-IR LFs. As 
AGNs are unimportant at log(Li5/L Q ) < 10.6, we extended 
the SF LFs to lower luminosities using the LFs from the 
MIPS sample. Evidently, the SF LFs were well described by 
Schechter profiles (Eq. [1]), similar to LFs at shorter wave- 
lengths (§ 1). We found the best-fit Schechter function param- 
eters with MPFIT and estimated their uncertainties with 1000 
Monte Carlo realizations (Table 4). 

We confirmed the strong, almost pure luminosity, evolu- 
tion in the SF LFs between < z < 0.7, as one can simply 
shift the local 15 fim LF by ^0.6 dex in luminosity to fit 
the total 15 fim LF at z ~ 0.7 (Fig. 6). Note that the local 
15 /tm LF of Xu (2000) exhibits a shallow bright-end slope, 
because the author did not differentiate between AGN and SF. 
On the other hand, optically selected AGNs were removed 
from the local (z < 0.3) 8 fim LF of Huang et al. (2007), and 
the authors employed a Schechter function to describe the 
SF LF. Quantitatively, the integrated 8 /im luminosity den- 
sity (Table 4) decreased by ^60% from z ~ 0.7 to z ~ 0.3. 
Using the SFR-L 8 relation of Huang et al. (2007) 15 , SFR = 
Ls/(L78 x 10 9 Lq) Mq yr -1 , this luminosity density gives 
a SFR density of 0.043 ±0.021 M yr" 1 Mpc" 3 at z ~ 0.7, 
where we have assumed a 50% calibration uncertainty (cf. Wu 
et al. 2005; Reddy et al. 2010). We defer a further estimate on 
the SFR density to § 4.3. 

Our AGN LFs agree with the obscuration-corrected AGN 
bolometric LF (Hopkins et al. 2007). The bolometric LFs 
were determined using observed AGN LFs from X-ray to 
mid-IR between < z < 6; obscured AGNs 16 were accounted 
for using the observed luminosity-dependent AGN absorp- 
tion distribution (Ueda et al. 2003) and assuming equal num- 

15 Converted to Chabrier (2003) IMF. 

16 Throughout the paper, we use "type-1" and "type-2" to describe op- 
tical spectral types and we use "obscured" and "unobscured" to distinguish 
whether or not the X-ray absorption column density is greater than 1 22 cirr 2 . 



bers of Compton-thick objects (Nn > 10 24 25 cm 2 ) and AGNs 
with A^h = 10 23-24 cm" 2 . The model including Compton-thick 
AGNs provides a better fit to the X-ray background spectrum. 
We ^-corrected the bolometric LF to rest-frame 8 /im and 
15 /im with the luminosity-dependent AGN SED templates 
of Hopkins et al. (2007). The ^-corrected LFs fit our AGN 
LFs at both wavelengths remarkably well (green solid curves 
in Fig. 6). The AGN contribution is less prominent in the 
8 /im LF than in the 15 /im LF because the PAHs increase 
the SF/AGN contrast. This result is reassuring, as mid-IR 
spectroscopy should provide the most complete AGN identi- 
fication. The obscuration correction factor for the bolometric 
LF is about 2.7 (i.e., obscured:unobscured ~ 1.7:1) in the lu- 
minosity range of the IRS sample. The IRS-selected AGNs, 
therefore, should be predominantly obscured sources. This is 
confirmed by (1) their low X-ray detection rate (§ 5.2), and 
(2) the rarity of broad emission lines in their optical spectra. 

Matute et al. (2006) built 15 /im LFs for type-1 AGNs at z 
= 0. 1 and 1 .2 and for type-2 AGNs at z = 0.05 and 0.35 using 
optically selected samples. We evolved their LFs to z = 0.7 
and combined the number densities of both optical types 17 
(blue dot-dashed curve in Fig. 6). The result agrees well 
with our AGN 15 /im LF. Assuming a flat AGN mid-IR SED 
in vLv, it also fits our AGN 8 fim LF. Even better matches 
were observed when we kept the emission powered by SF in 
the AGN/SF composite systems, imitating the situation faced 
by Matute et al., who were unable to decompose the mid-IR 
emission. Nevertheless, since evolving the LFs involved large 
extrapolation for the type-2 AGN LF, which dominated the 
AGN number densities, this should be regarded as only tenta- 
tive evidence that optically selected AGN samples are as com- 
plete as mid-IR selected AGN samples at \og(L\slLo) > 10.5. 

The comoving space density of AGNs rises above that of 
star-forming galaxies at log(Li5/L Q ) > 10.8, and it is dom- 
inated by type-2 AGNs at log(Li 5 /L Q ) < 11.5. Only six of 
the 20 IRS-selected AGNs with optical spectra show broad 
emission lines; the rest only show narrow emission lines. The 
predominance of type-2 AGNs over star-forming galaxies at 
high luminosities is supported by the stacked optical spec- 

17 Since all of their evolution models produced similar results in the bright 
end, we show only the result from the models with evolving faint-end slopes. 
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trum of 24 /mi sources at Lig ~ vL v ^ > 10" L Q between 
0.4 < z < 0.7, as presented in Caputi et al. (2008). The spec- 
trum shows narrow emission lines and an [O lll]5007/H/3 ra- 
tio of ~3 ; indicating AGN photoionization. In contrast, the 
emission-line ratios of 24 /im sources at lower luminosities 
are more consistent with star-forming regions. 

4.3. 8—1000 fjm Luminosity Function 

Since the bolometric output of obscured SF emerges mostly 
as cold dust emission (T ~ 40 K), the far-IR wavelengths most 
directly measure the energetics of active SF, and subsequently, 
the SFR. Another advantage of the far-IR regime is that the 
AGN contribution should be minimal, since the average AGN 
SED drops sharply beyond ^20 /im. To take a glimpse of 
the SFR function, we estimated the 8-1000 /im LF at z ~ 0.7 
using a 70 /im selected sample. 

Again, we used the 1/V max estimator. Our sample consisted 
of 70 /mi-detected sources above 7 mJy between 0.6 < z < 
0.8. The flux threshold is the 3er detection limit for the nomi- 
nal coverage of 100 for the S-COSMOS field (for a coverage 
map see Frayer et al. 2009). The sample area is 1.66 deg 2 , the 
same as the 24 /im sample. 

Symeonidis et al. (2008) presented a simple equation to 
convert observed luminosities at 24, 70, and 160 /mi to Lir. 
This equation 18 produces similar results as more sophisticated 
SED modeling for a wide redshift range (a ~ 0.06 dex; Kar- 
taltepe et al. 2010). We thus computed Lir using this equation. 
For the 70 /im sources undetected at 160 /im, we assumed an 
average L160/L70 ratio 19 of 4, according to the S-COSMOS 
median stacking analysis (Lee et al. 2010). We adopted the 
completeness curve for coverages greater than the nominal 
value of 100 to correct for source extraction incompleteness 
below L70 = 16 mJy (Frayer et al. 2009). The W max LF is 
shown in Figure 7 and is tabulated in Table 5. 

We compared our LF with that of Magnelli et al. (2009) 20 , 
who also used 70 /im data to estimate the bright end of the 
8-1000 /im LF. We interpolated their LFs to z = 0.7. De- 
spite an offset, the bright-end slopes are consistent down to 
the lowest comoving density accessible to their survey (their 
largest field is 0.25 deg 2 ). Note that Magnelli et al. fixed the 
bright-end slope to that of the local 8-1000 /im LF (<j> oc L~ 2,2 ; 
Sanders et al. 2003). The ^0.2 dex offset is due to the differ- 
ence in deriving luminosities. The SED templates of Chary 
& Elbaz (2001), which Magnelli et al. used to compute Lir, 
gives /^-corrections (LirAvL^, ^ b "') 20-40% lower than those 
from the equation of Symeonidis et al. (2008). 

Magnelli et al. (2009) showed that the 8-1000 /im LFs are 
consistent with the 15 /im LFs once the latter are ^-corrected 
with their stacking analysis results. This consistency is ex- 
pected because both mid-IR and total IR probes obscured 
SF. This result, however, is moot because their 15 /im LFs 
are heavily contaminated by AGNs but their far-IR-derived 
8-1000 /im LFs are probably not. Since our SF 15 fim LF 
has a steeper bright-end slope than the total 15 /an LF, we 
re-examined the relation between mid-IR and total IR LFs. 

Previous stacking analyses and various luminosity- 
dependent IR SED libraries have consistently shown that the 

18 The (1+z) term must be dropped if one does not de-redshift the observed 
flux densities. 

19 Interpolating the stacked F160/F 70 values according to F24 or F70 
fluxes makes almost no difference in terms of the overall shape of the LF. 

20 For comparisons with other determinations of the LF, refer to Fig. 12 of 
their paper. 



extrapolation from L15 to Lir is independent of L15 or Lir 
(e.g., Zheng et al. 2007b; Lee et al. 2010; see, e.g., Fig. 8 of 
Le Floc'h et al. 2005). Using the average Lir/Ljs ratio from 
the S-COSMOS median stacking analysis (<Lt R /Li 5 > = 19 at 
0.6 < z < 0.8; Fig. 8 Lee et al. 2010) 21 , we converted the best- 
fit SF 15 /tm LF to an 8-1000 /mi LF (black dashed curve in 
Fig. 7). For consistency, we calculated the luminosities of 
the stacked galaxies with the equation of Symeonidis et al. 
(2008), although they agreed with those from the SED fitting 
of Lee et al. (2010). As we suspected from the steep bright- 
end slope of the 15 /im LF, this /T-corrected LF significantly 
underestimates the 8-1000 /im LF. 

In the median stacking analysis that we used, the average 
bolometric conversion factor has included both detections and 
non-detections at 70/160 /im. Hence, the observed mismatch 
is not due to an underestimated Lir/Lis. In fact, previous 
mean stacking analyses gave much smaller conversion factors 
(<Lir/Li 5 >~ 8; Zheng et al. 2007b; Magnelli et al. 2009), 
because 70/160 /mi detections were excluded. 

An important factor that we have neglected so far is the 
scatter in the L15-L1R correlation. We show the dispersion 
in Lir/Lli by plotting F70/F24 against F24 (Fig. 8). For 
objects selected within a narrow redshift bin, like ours, the 
effect of differential /T-correction is negligible. Since F70 
strongly correlates with Lir and L24 is approximately 15 /im 
flux at the rest-frame, L70/L24 is equivalent to Lir/Lis. The 
70 /im detection fraction increases with 24 /im flux, with a 
profile similar to an error function (bottom panel of Fig. 8) 22 . 
Given that Lee et al.'s stacking analysis results represent the 
average L70/L24 for all objects within a given 24 /im flux 
bin, the increase of 70 fim detection fraction with 24 /im flux 
can be understood if the L70 / L24 distribution is Log-Normal 
with a = 0.23 dex. Admittedly, this measurement is only valid 
above L24 > 0.3 mJy, because of the small detection fraction 
at lower fluxes. But our conclusion is insensitive to the scatter 
at low 24 /im fluxes. 

The dispersion is mostly intrinsic and reflects the variations 
in IR SEDs of SF. The measurement errors (^0. 1 dex), which 
include both random errors in the photometry and the sys- 
tematic calibration uncertainties at both 24 and 70 /im (Engel- 
bracht et al. 2007; Gordon et al. 2007), are low compared to 
this dispersion. From the stacking analysis results at adjacent 
redshift bins, we estimated that the spread due to our finite 
redshift bin is about ±0.05 dex. After subtracting quadrati- 
cally the measurement uncertainties and the spread due to the 
redshift bin, we found that the intrinsic dispersion of the flux 
density ratio is 0.2 dex. Note that nearby IR galaxies show 
the same amount of dispersion in the Lir/L^ ratios (Chary & 
Elbaz 2001). 

To evaluate the impact of this dispersion to the LF transfor- 
mation from 15 /im to 8-1000 /tm, we generated a mock cat- 
alog of galaxies that follows the Schechter profile of the SF 
15 /im LF and the redshift distribution of the actual catalog. 
We then applied inverse ^-corrections to obtain their 24 fim 
fluxes, and we randomly assigned their 70 /im fluxes using 
the stacking analysis results and the dispersion of L70/L24. 
The 160 /im fluxes were fixed to the average L160/L70 ratio 
from stacking. Again we calculated Lir using the equation 

21 We computed Lir for each 24 fim flux bin using the stacking results and 
the equation of Symeonidis et al. (2008), which gave consistent results as the 
full SED modeling of Lee et al. (2010); we computed L\s by /^-correcting the 
24 /im fluxes in the same way as in § 4. 1 . 

22 As AGNs dominate above F24 > 1 mJy, we used the AGN-subtracted 
IRS sample to calculate the 70 fim detection fraction at F24 > 0.7 mJy. 
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FIG. 6. — (a) 15 /im and (fo) IRAC 8 /tm LFs at z ~ 0.7. The IRS decomposed AGN and SF LFs are shown as blue filled circles and red open circles, respectively. 
The SF LFs were extended to lower luminosities with the MIPS sample down to F24 = 60 fiiy (black/grey squares show data points below/above the luminosity 
limit of the IRS sample). The best-fit Schechter functions of the SF LFs are shown as black solid curves, with the shaded areas delimiting the lcr uncertainties 
from Monte Carlo realizations. As a comparison with previous studies, in the left panel we show the 15 /im LF at 0.6 < z < 0.8 of Le Floc'h et al. (2005) 
{grey crosses with error bars). The AGN LFs match the 15 fim LF of optically selected type-1 and type-2 AGNs (blue dashed curve; Matute et al. 2006) and the 
^-corrected obscuration-corrected AGN bolometric LFs (green solid curves', Hopkins et al. 2007). To show LF evolution, we included the local 15 fim and 8 /im 
LFs (purple pentagons and dotted curves; Xu 2000; Huang et al. 2007). The top abscissa indicates corresponding 24 fim fluxes at z = 0.7. The knees (L*) of the 
LFs are marked by stars. 
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FIG. 7.— 8-1000 fim LF at z ~ 0.7. The 1/V ma x data points from the 70 fim 
selected sample are black squares and grey circles (before and after remov- 
ing power-law AGNs, respectively). The blue lines extend to the 1/V ma x data 
points before correcting for catalog incompleteness to illustrate the uncertain- 
ties below 10 12 L Q . The green dashed curve is the best-fit SF 15 /im LF, bolo- 
metrically corrected adopting the mean Lir/L^ ratio of 19 from the stacking 
results (Lee et al. 2010). The solid curve shows the synthetic LF based on 
the same SF 15 /im LF but incorporated scatter in bolometric correction (see 
§ 4.3 for details). For comparisons, we also included the 8-1000 fim LF of 
Magnelli et al. (2009) (red dotted curve) and the local 8-1000 fim LF (purple 
pentagons; Sanders et al. 2003). The knees (L*) of the LFs are marked by 
stars. 

of Symeonidis et al. (2008). As a convolution of a Gaussian 
and a Schechter function, the synthetic LF is best described 
by a modified Schechter profile (Eq. [2]), although Schechter 
function can also provide an acceptable fit. The best-fit model 
is shown as the black solid curve in Fig. 7 and its parameters 
are listed in Table 4. With a shallower bright-end slope, it is 
now consistent with the 70 /jm-derived LF and its extrapo- 
lation to 10 13 Lq agrees with that of sub-milimeter galaxies 
(Chapman et al. 2005). 
The synthetic LF also extends the 8-1000 fim LF below the 



"knee". Integrating the synthetic LF we obtained a total bolo- 
metric luminosity density of 7 x 10 8 L Q Mpc~ 3 . Using the 
Kennicutt (1998) calibration, SFR = L IR /(10 10 L ) M yr-\ 
we obtained a SFR density of 0.07 ± 0.02 M Q yr" 1 Mpc -3 , 
where the error is dominated by the ^30% calibration uncer- 
tainties. This result agrees well with previous estimates using 
non-IR tracers (see the compilation of Hopkins & Beacom 
2006). LIRGs/ULIRGs contribute about 54%/7% of the total 
SFR density at z ~ 0.7. The ULIRG contribution to the total 
SFR density would have been underestimated by a factor of 
four if we had simply converted the SF 15 /im LF using the 
mean Lir/L^. 

In summary, only after incorporating both the average 
Lir/Lis and its dispersion were we able to correctly predict 
the 8-1000 fim LF from the SF 15 fim LF. Note that our 
result is insensitive to the shape of the Lir/Lis distribution: 
whether or not it is Log-Normal, this mismatch in the bright- 
end slope can be reconciled as long as ^16% of the sources 
with L15 >2x 10 10 Lq have Lir/Lis ratios 1.6 times higher 
than the mean value. We suggest two possible physical ori- 
gins of these galaxies: 

1 . These galaxies show the hottest IR SEDs because they 
are extreme starbursts. It is known that star-forming 
galaxies form a SFR-M* sequence with an intrinsic dis- 
persion of <0.3 dex (Noeske et al. 2007; Zheng et al. 
2007a; Elbaz et al. 2007; Daddi et al. 2007). Since the 
SFRs in these studies were derived mostly using 24 fim 
data, the dispersion in the true SFR-M* sequence must 
be larger than the current estimates given the ~0.2 dex 
dispersion in the ratio of 70 fim and 24 /im-derived 
SFRs. In analogy, this dispersion can explain why the 
SFR Function has a much shallower bright-end slope 
than the stellar Mass Function of star-forming galax- 
ies (Bell et al. 2007; Ilbert et al. 2010; Boissier et al. 
2010). The dispersion could naturally arise if there is 
a large range of SF efficiency (SFR/Mm) (Scoville & 
Soifer 1991; Daddi et al. 2010) and/or molecular gas 
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F24 [mJy] 

FIG. 8. — 70/24 fim flux density ratio vs. 24 fim flux for galaxies at 
0.6 < z < 0.8. Individual detections above 7 mJy are shown as grey filled 
circles. The green data points highlight the 30 IRS sources detected at 70 (im; 
and the big circles show AGN-corrected results for the 19 IRS sources with 
F24(SF) > 0.7 mJy, color-coded according to the AGN contribution (pagn; 
see the legend). The grey dashed line indicates F70 = 7 mJy, the 3<r detection 
limit. The stacking analysis results of Lee et al. (2010) are shown as black 
diamonds connected by a solid curve; and the error bars show the expected 
variation due to ^-correction within 0.6 < z < 0.8. In the bottom panel, 
the solid histogram shows the fraction of 24 fim sources detected at 70 fim. 
To show the effect of the catalog incompleteness correction, the dotted his- 
tograms show the detection fraction without correcting for this incomplete- 
ness. AGN corrections based on IRS spectra were applied for bins above 
F24 > 0.7 mJy. The blue curve shows the expected fraction of 70 fim de- 
tected sources if the F10/F2A ratio is a Log-Normal distribution with <j = 
0.23 dex and means centered on the stacking results. The dashed curves in 
the upper panel mark this dispersion (±1<t). 

fraction (Mui/M*) in star-forming galaxies at a given 
stellar mass. 

2. The rest-frame ~40 /im emission in these galaxies 
could be powered by AGNs, similar to Mrk231 
(Veilleux et al. 2009; van der Werf et al. 2010). In 
fact, if we remove power-law AGNs (a\ > -0.2; § 5.2) 
from the 70 /im selected sample, then the 8-1000 LF 
becomes more consistent with the SF 15 /im LF as it 
is significantly reduced at Lvr > 10 12 Lq (grey circles 
in Fig. 7). Although these data points should be con- 
sidered as lower limits because the IR luminosities in 
these AGNs could also be powered by obscured SF, 
it is intriguing that almost all of the IRS sources with 
F1Q/F2A values 0.23 dex above the stacking results 
contain AGNs (Fig. 8). 

5. DISCUSSION 

5.1. Co-Evolution ofBHs and Galaxies in LIRGs 

Previous multi-wavelength estimates of the AGN and SF 
LFs have shown that the volume-averaged BH accretion his- 
tory closely follows the volume-averaged SFR history at z < 2 
(e.g., Shankar et al. 2009; Zheng et al. 2009). The ratio of 
the two, M B h(z)/SFR(z) ~ 5-8 x 10" 4 , matches not only the 
present-day ratio of the BH mass density and stellar mass 



density (psn/Pstar — 1.5 x 10~ 3 ) but also the normalization 
of the local M B H-M b uige relation (M B H/M bu i ge ~ 1 .4 x 10" 3 at 
Mbuige = 5 x 10 10 M Q ; Haring & Rix 2004), provided that part 
of the stellar mass is recycled into the interstellar medium. 
These observations strongly suggest a co-evolution of galax- 
ies and BHs, but are they co-evolving in the same objects and 
at the same time? 

With the AGN/SF spectral decomposition results (§ 3), we 
can directly address this problem using the 48 galaxies with 
L 15 > L 1 ™' = 3 x 10 10 L & and between 0.6 < z < 0.8. The 
15 /.mi luminosity limit corresponds to L lR ~ 6 x 10" Lq 
or SFR ~ 60 Mq yr" 1 for star-forming galaxies and L b oi — 
4 x 10 11 L Q or M B h = 0.24// .i Mq yr" 1 for AGNs, where 
/o.i = e /(l ~ e ) f° r £ = 0.1, which is the BH radiation efficiency. 
Within the sample, 31 sources have SF luminosities (iff) 
greater than Lj™ 1 , 24 have AGN luminosities (L^ 5 GN ) greater 
than L 1 ™ 1 , and 7 have both L?f and Lf 5 GN above L 1 ^. 

Intense BH accretion and SF therefore coexist in ^23% 
(7/31) of the (U)LIRGs (SFR > 60 M Q yr" 1 ). Assuming 
every (U)LIRG experiences such an AGN phase in its life- 
time, the AGN duty cycle is 0.23 ± 0.09 23 for (U)LIRGs 
at z ~ 0.7. If we assume the observed luminosity distribu- 
tion evenly samples the light curve of a typical star-forming 
galaxy or an AGN, then the mass ratio of BH growth and SF 
can be obtained by integrating luminosities for each compo- 
nent over the 31 galaxies with SFR > 60 Mq yr" 1 . The re- 
sult is AM BH /AM slar ~ 1.7 x 10~ 3 //o.i. For the 7 compos- 
ite galaxies with both components above the luminosity limit, 

we measured ][]Mbh/E sfr - 4 - 6 x 10 _3 //o.i. Multiply- 
ing by the AGN duty cycle (7/31), the time-averaged ratio is 
Mbh/SFR~ 1.0 x 10- 3 //o.i. 

Our estimates are consistent with the local MfiH-A^buige 
relation (Haring & Rix 2004), supporting an invariable 
MfiH-ZWhost relation since z ~ 1 (Jahnke et al. 2009; Bennert 
et al. 2010). For (U)LIRGs, this result also rules out co- 
evolution scenarios in which the bulges and BHs grow in the 
same objects and at the same time or SF and BH accretion 
occur in different events. Instead, it favors a co-evolution 
with a time offset between SF and BH accretion, in which 
BH growth and SF happen in the same event but the former 
has a much shorter lifetime than the latter. 

For the 17 AGNs with Lf GN > L'f but L$ < L 1 ™', we 

measured ^Mbh/X^FR — 1.2 x 10 -2 / /o.i, almost an or- 
der of magnitude higher than the normalization of the lo- 
cal MfiH-A^buige relations. However, since the galaxy co- 
moving density increases dramatically as SFR decreases, the 
AGN duty cycle in these galaxies is much lower than that 
of the (U)LIRGs with SFR > 60 M yr" 1 . With zCOS- 
MOS spectra, Silverman et al. (2009) found that on average 
Mbh/SFR ~1.9x 10" 2 for X-ray selected AGNs at z < 1. 
These AGN hosts have a mean SFR ~ 10 Mq yr -1 . The 
four times higher Mbh/SFR than what we measured in the 
composite (U)LIRGs also suggests a smaller AGN duty cy- 
cle, probably as a consequence of the longer lifetime for the 
SF. 

5.2. Comparison with Other AGN Identification Methods 

Mid-IR spectral analysis offers arguably the best method to 
identify AGNs. In most cases, however, one has to rely on 

23 Poisson error. 
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FIG. 9. — Our mid-IR spectroscopical selection of AGNs compared with major AGN identification methods. Symbols are the same as in Fig. 4. The median 
errors are shown at the top left corner in each panel, (a) IR excess (q) vs. radio luminosity. Downward arrows mark the radio sources undetected at 70 /im. The 
grey dotted and dashed lines indicate the radio selection thresholds (q < 1.64 or Li 4GHz > 6.3 X 10 23 W Hz~' ; Yun et al. 2001; Hickox et al. 2009). The red 
dot-dashed line shows the IR luminosity limit of the 70 /im sources at z ~ 0.7 (Lir > 5 X 10" Z.q). (b) Optical— mid-IR color vs. L15. X-ray AGNs (stars) show 
bluer colors (i.e., less obscured) than non-X-ray AGNs. Sacchi et al. (2009)'s Type-l/Type-2 AGN separation line is the dashed line (F24/F(r) = 100). The three 
highly obscured X-ray AGNs are labeled (§ 5.3). (c) IRAC color-color diagram. The dashed lines delimit the AGN locus of Stern et al. (2005). (d) Rest-frame 



near-IR spectral index (a^ = \og(uLj"" si / uL^" st )) vs. rest-frame 5 fim excess in magnitude. The 5 (im excess is the magnitude difference between the observed 
8 fim flux and the stellar blackbody emission. Most AGNs would be identified if one selects objects with either o& > -0.2 or 5 /^m excess < — 1 .4 mag. 



photometric data to select AGNs, since it is difficult to carry 
out mid-IR spectroscopy for a large number of sources. We 
compared our mid-IR decomposition results with commonly 
used photometry-based AGN identification methods to select 
the most efficient methods for future studies. We emphasize 
that since our comparisons are limited to the IRS sample, i.e., 
24 /m sources above 0.7 mJy and between 0.6 < z < 0.8, we 
cannot comment on the contaminations by sources below the 
24 /im flux limit and/or at other redshifts. 

From our spectral decomposition of the 53 IRS sources 24 , 
we estimated both the AGN luminosity and its contribution to 
the total mid-IR SED. For 31 sources, a power-law AGN was 
required to fit the mid-IR spectrum. In four such sources, the 
AGN luminosity (Lff N ) fell below 3 x 10 10 L (Lj™'). Al- 
though these four sources are possibly real AGNs, we ignored 
them in this discussion because the sample is incomplete be- 
low the luminosity limit. We thus have 27 "AGN hosts" with 
Lff N >L l {f and 21 'AGN-dominated sources" with p AGN > 

24 These include the 5 objects outside of the redshift bin, 0.6 < z < 0.8. 



50%. All of the AGN-dominated sources are also AGN hosts. 

41 IRS sources have radio counterparts from the VLA- 
COSMOS survey (Schinnerer et al. 2007). Three can be 
classified as radio AGNs (Fig. 9a), as their rest-frame ra- 
dio luminosities are greater than 6.3 x 10 23 W Hz" 1 (e.g., 
Hickox et al. 2009). An alternative selection method is based 
on the IR-radio correlation (Helou et al. 1985). Yun et al. 
(2001) suggested that sources with q parameters five times 
less than the mean value (i.e., q < 1 .64) are AGNs, where q = 
log(LrR/L )-log(Li.4GHz/lO 14 W Hz" 1 ). Our most radio lumi- 
nous source is also the only radio-excess object with q = 1 .54 
(see Sargent et al. 2010). We calculated Lir in the same way 
as in § 4.3. The 11 radio sources without 70 /im detection 
are marked by downward arrows. Radio power/radio excess 
identified only 1 1/4% (3/27,1/27) of the AGN hosts. 

14 sources have X-ray counterparts from the XMM- 
COSMOS survey (stars in Fig. 9; Hasinger et al. 2007; Brusa 
et al. 2010), including the 3 radio-selected AGNs. These are 
bona fide AGNs, since their rest-frame X-ray luminosities 
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(Lx = £o.5-iokev) are all greater than 10 43 erg s 1 given a pho- 
ton index of V = 1 .7. The X-ray selection, however, identified 
only 52% (14/27) of the AGN hosts and 62% (13/21) of the 
AGN-dominated sources. The central 0.9 deg 2 of the COS- 
MOS field was covered by the much deeper Chandra COS- 
MOS survey (Elvis et al. 2009, Civano et al. 2010, in prep.). 
The Chandra detection limit was /o.s-iokev = 5.7 x 10~ 16 erg 
s" 1 cm" 2 , corresponding to Lx = 10 42 erg s" 1 at z ~ 0.7. Eight 
XMM-selected AGNs were inside the Chandra coverage, and 
they were all in the Chandra catalog with consistent fluxes. 
Additionally, Chandra detected two pure star-forming galax- 
ies (/?agn = 0) with Lx < 1.5 x 10 42 erg s" 1 that were not in 
the XMM catalog. We attributed the X-ray emission of these 
two objects to SF activity. Although 11 of the 13 non-X-ray 
AGNs were observed by Chandra, none was detected. 

The Spitzer/IRAC -based AGN diagnostics select sources 
with power-law SEDs, presumably dominated by AGN heated 
hot dust emission. The IRAC color-color selection of Stern 
et al. (2005) identified 67% (18/27) of the AGN hosts and 
71% (15/21) of the AGN-dominated sources (Fig. 9c). The 
near-IR spectral index between rest-frame 2 and 5 /im, a\ = 
\og{L^/L2), is even more efficient. Since the effective wave- 
lengths of the IRAC 3.6 and 8.0 /tm filters correspond to rest- 
frame 2. 1 and 4.5 /im at z ~ 0.7, we computed a\ using [3.6] 
and [8.0] magnitudes directly. Applying a\ > -0.2, we iden- 
tified 89% (24/27) of the AGN hosts and 100% (21/21) of the 
AGN-dominated sources. The rest-frame 5 /im excess selec- 
tion works equally well, as indicated by its tight correlation 
with a, (Fig. 9d). We computed the excess by subtracting 
the extrapolated stellar blackbody emission from the observed 
IRAC 8 /im magnitude. The extrapolation used the best-fit 
Bruzual & Chariot (2003) model to the FUV-to-K s SED (§ 3). 
The 2-to-5 /im spectral index and the rest-frame 5 /im excess 
are, therefore, the best photometric AGN diagnostics, as their 
results are most consistent with the mid-IR spectral decompo- 
sition results. 

5.3. Compton Thick AGNs? 

The remarkable agreement between mid-IR LFs of IRS- 
selected AGNs and obscuration-corrected AGN bolometric 
LFs implies that we may have achieved a complete AGN 
identification (§ 4.2; Fig. 6). We have confidence in this 
result because of our high spectroscopic redshift complete- 
ness (79%; § 2.1) for X-ray-detected bright 24 pm sources 
(F24 > 0.7 mJy). Since most unobscured AGNs with Lf° N > 
Lj™' should be detected in X-ray at z ~ 0.7 and most ob- 
scured AGNs should not have strange optical SEDs, the red- 
shifts of the AGNs in the parent sample should be accurate. 
Recall that Compton-thick AGNs were counted in the bolo- 
metric LFs, using the A^h distribution of Ueda et al. (2003) 
and assuming as many Compton-thick objects as AGNs with 
N H = 10 23 " 24 cm" 2 (Hopkins et al. 2007). Following the same 
A^h distribution, we estimated that roughly 35%, 50%, and 
15% of the IRS-selected AGNs were unobscured, Compton- 
thin (A^h = 10 22 " 24 cm" 2 ), and Compton-thick, respectively. 

The optical-to-24/im color has been widely used to select 
obscured AGNs (Fiore et al. 2008, 2009; Dey et al. 2008). Our 
non-X-ray AGNs appear more obscured than X-ray AGNs as 
they show redder r AB - [24] ab colors; and the two groups are 
roughly separated at F24/F(r) = 100, the type-l/type-2 divid- 
ing line suggested by Sacchi et al. (2009) (Fig 9b). Since six 
of the 1 1 X-ray AGNs and zero of the nine non-X-ray AGNs 



with optical spectra show broad emission lines, we could fur- 
ther assume that most of the unobscured AGNs were detected 
in X-ray. Hence, only ^23% of the obscured AGNs are de- 
tected in X-ray. Are there any Compton-thick ones in these 
X-ray AGNs? The three X-ray AGNs with the highest 15 /im 
luminosities (MIPS 42682, 45149, & 51645; see the last three 
panels in Fig. 3) show the reddest F24/F{r) colors as well as 
the highest X-ray hardness ratios, HR = (H-S)/(H+S) > 0.5 
(others have HR < 0), where H and S are the counts in 
the 2-10 and 0.5-2 keV energy bands, respectively. Two 
of the three have optical spectra, and both show only nar- 
row emission lines. Although unlikely Compton-thick, these 
sources are heavily obscured quasars, as HRs above 0.5 indi- 
cate A^h > 10 23 cm" 2 and their absorption-corrected X-ray lu- 
minosities are much greater than 10 44 erg s" 1 (Mainieri et al. 
2007). In fact, MIPS 45149 (XID = 70, HR = 0.7) was in- 
cluded in Mainieri et al.'s X-ray spectral analysis; they mod- 
eled its X-ray spectrum using an absorbed power-law with 
A^h = 1 .7 x 10 23 cm" 2 and an intrinsic rest-frame luminosity of 
4 n ?-i0kev = 6.8 x 10 44 erg s" 1 . In the IRS sample, MIPS 45 149 
is also the most heavily extinct X-ray AGN in the mid-IR 
(r 9 A 7 GN = 2.1). 

It is natural to suspect that the four objects with the high- 
est mid-IR extinctions for the AGN component (MIPS 30035, 
48730, 49150, 51523), t 9 a ™ = 3.4-5, are the Compton-thick 
sources (see Fig. 3), since (1) all of the four were undetected 
in X-ray, (2) 80% (21/27) of the IRS-selected AGNs have 
t^ n < 1, and (3) the highest t^ n measured in the X-ray 
AGNs is only 2.1. Note, however, that these optical depths 
correspond to only A^h = 1-1 .6 x 10 23 cm" 2 if one assumes the 
Galactic dust-to-gas ratio and the extinction law of Rieke & 
Lebofsky (1985). By definition, without X-ray data on these 
objects, we are unable to classify these sources as Compton- 
thick AGNs. 

6. CONCLUSIONS 

The mass assembly history of galaxies and their BHs are 
imprinted in time-resolved IR LFs. Our understanding of both 
processes have been hindered by the degeneracy of BH ac- 
cretion and obscured SF in the LFs, especially in the mid- 
IR regime. In this paper we used Spitzer mid-IR spectra to 
break this degeneracy and decompose mid-IR LFs between 
SF and AGNs at z ~ 0.7. Further, we built a 8-1000 pm LF 
with 70 /im data to evaluate the high-end of the SFR func- 
tion and to examine the relation between mid-IR LFs and 
the SFR Function. Finally, we discussed the implications of 
our analysis to BH-galaxy co-evolution scenarios, AGN se- 
lection methods, and the elusive Compton-thick AGN popu- 
lation. Our conclusions can be summarized as follows: 

1. The SF mid-IR LFs have steep bright-end slopes. Their 
profiles are well described by the classic Schechter 
function over three decades in comoving density, simi- 
lar to LFs at shorter wavelengths. 

2. The AGN mid-IR LFs have shallower bright-end 
slopes. AGNs are thus responsible for the bright-end 
excess of mid-IR LFs. The AGN mid-IR LFs, once K- 
corrected, match the obscuration-corrected AGN bolo- 
metric LF, which was determined with the observed 
AGN multi-wavelength LFs after accounting for ob- 
scured AGNs (both Compton-thin and Compton-thick) 
with the AGN absorption distribution. The Nh distri- 
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bution predicts that 15% (—4) of the AGNs in our IRS 
sample are Compton-thick. 

3. Although our AGNs are mid-IR selected, the LFs are 
also consistent with extrapolations from previous de- 
terminations with optically-selected AGNs. This ten- 
tatively implies that optical spectroscopy can identify 
most of the obscured AGNs, at least for those with high 
IR luminosities. 

4. The 8-1000 /im LF derived from a 70 pm selected 
sample shows a shallower bright-end slope than the 
bolometrically-corrected SF 15 /im LF. The mismatch 
can be reconciled only after incorporating the intrin- 
sic dispersion in the Z45-Z4R correlation, especially 
the ^16% of 24 /im sources with high L\r/L\5 ratios 
(> 1.6x higher than the mean). These sources could 
be either extreme starbursts or SF/AGN composite sys- 
tems in which bulk of the far-IR emission is powered 
by BH accretion. We further suggest that the disper- 
sion in the SFR-M* sequence explains why the SFR 
Function has a much shallower bright-end-slope than 
the stellar Mass Function of star-forming galaxies and 
the dispersion naturally arises if there is a large range of 
SF efficiency (SFR/MH2) and/or molecular gas fraction 
in star-forming galaxies at a given stellar mass. 

5. Intense BH accretion accompanies about a quarter of 
a LIRG's lifetime when the SFR soars above 60 M© 
yr" 1 . Throughout the LIRG's lifetime, the growth in BH 
mass is ~0.1% of the mass of newly formed stars, con- 
sistent with the local MBH-^buige relation, provided that 
part of the stellar mass is recycled into the interstellar 
medium. These results support a constant MBH-^host 
relation since z <~ 1 and favor co-evolution scenarios 
in which BH growth and SF are triggered in the same 
event but the former spans a much shorter lifetime than 
the latter. 

6. X-ray selection (Lx > 10 42 erg s" 1 ) misses half of 



the mid-IR identified AGNs. The majority of the 
X-ray AGNs are unobscured, despite that the three 
most mid-IR-luminous X-ray AGNs appear highly ob- 
scured (A^h > 10 23 cm" 2 ). The majority of the non- 
X-ray AGNs are obscured. We suspect that the four 
sources with the highest silicate extinctions, as mea- 
sured in their power-law components (r^ N = 3.4-5), 
are Compton-thick AGNs. All of the four are unde- 
tected in X-ray. 



7. Among the photometric AGN identification methods, 
near-IR spectral index and IRAC excess are most con- 
sistent with our mid-IR spectral decomposition analy- 



Although our results are limited to a narrow redshift bin 
(0.6 < z < 0.8), the above conclusions could be extended to 
other redshifts. The evolution of IR LFs at < z < 3 based on 
Spitzer and Herschel data will be the subject of a future paper 
(Le Floc'h et al. 2010, in prep.). Extending the AGN/SF 
decomposition of IR LFs to lower luminosities and for a wider 
redshift range will greatly improve the measurements of the 
growth rates of galaxies and their BHs. Together with better 
constraints on the redshift and mass-dependent normalization 
in the mass-growth-rate relations (i.e., the specific SFRs and 
the BH Eddington ratios), we will finally be able to retrace the 
evolutionary paths of galaxies and their BHs and to witness 
how the Magorrian relation emerged. 
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Decomposing Star Formation and AGN 



TABLE 1 
Spitzer IRS SAMPLE 



MIPS ID 
(1) 



Opt. RA 

(2) 



Opt. Dec 

(3) 



(4) 



zflag 
(5) 



SL1 

(6) 



LL2 

(V) 



LL1 

(8) 



F24 

(9) 



Z45 
(10) 



(ID 



PAGN 
(12) 



2176 


149.91600 


1.879972 


0.834 


spec-z 


20 


x 2 x 60 


15 


x 2 x 


120 


10 


X 


2 


X 


120 


5264 


149.80170 


2.384056 


0.704 


spec-z 


13 


x 2 x 60 


11 


x 2 x 


120 


7 


x 


2 


X 


120 


15123 


150.05206 


2.126712 


0.665 


spec-z 


10 


x 2 x 60 


10 


x 2 x 


120 


6 


X 


2 


X 


120 


16913 


149.87215 


2.289672 


0.700 


spec-z 


12 


x 2 x 60 


13 


x 2 x 


120 


8 


x 


2 


X 


120 


18902* 


149.75092 


2.469942 


0.657 


spec-z 


5 


x 2 x 60 


9 


x 2 x 


120 


7 


x 


2 


X 


120 


18984 


150.15953 


2.474338 


0.691 


spec-z 


15 


x 2 x 60 


14 


x 2 x 


120 


7 


x 


2 


X 


120 


20765 


150.04402 


2.643906 


0.697 


spec-z 


8 


x 2 x 60 


7 


x 2 x 


120 


5 


X 


2 


X 


120 


23930 


149.91074 


1.662555 


0.633 


spec-z 


8 


x 2 x 60 


12 


x 2 x 


120 


8 


x 


2 


X 


120 


26741 


149.71869 


1.849873 


0.677 


spec-z 


12 


x 2 x 60 


12 


x 2 x 


120 


6 


X 


2 


X 


120 


30035 


150.21017 


2.311710 


0.747 


spec-z 


8 


x 2 x 60 


10 


x 2 x 


120 


6 


X 


2 


X 


120 


31181* 


150.63902 


2.464456 


0.800 


spec-z 


4 


x2x60 


7 


x 2 x 


120 


6 


X 


2 


X 


120 


31635 


150.26385 


2.541551 


0.712 


spec-z 


11 


x 2 x 60 


11 


x 2 x 


120 


6 


X 


2 


X 


120 


36942 


150.42406 


2.123016 


0.699 


spec-z 


15 


x 2 x 60 


14 


x 2 x 


120 


7 


X 


2 


X 


120 


38331 


149.60471 


2.452618 


0.674 


spec-z 


8 


x 2 x 60 


10 


x 2 x 


120 


6 


X 


2 


X 


120 


38754 


149.73370 


2.575010 


0.706 


spec-z 


6 


x 2 x 60 


9 


x 2 x 


120 


5 


X 


2 


X 


120 


40291 


150.10332 


1.661929 


0.834 


spec-z 


17 


x 2 x 60 


14 


x 2 x 


120 


7 


X 


2 


X 


120 


40410* 


150.24252 


1.766390 


0.623 


spec-z 


4 


x 2 x 60 


8 


x 2 x 


120 


6 


X 


2 


X 


120 


42682* 


149.58641 


1.769320 


0.787 


spec-z 


2 


x 2 x 60 


2 


x 2 x 


120 


2 


X 


2 


X 


120 


42997 


150.66016 


1.863509 


0.624 


spec-z 


25 


x 2 x 60 


14 


x 2 x 


120 


9 


X 


2 


X 


120 


43504 


149.99924 


2.005987 


0.760 


spec-z 


3 


x 2 x 60 


4 


x 2 x 


120 


4 


X 


2 


X 


120 


44180 


150.21385 


2.188338 


0.799 


spec-z 


15 


x 2 x 60 


15 


x 2 x 


120 


8 


X 


2 


X 


120 


44883 


149.94170 


2.395768 


0.759 


spec-z 


10 


x 2 x 60 


12 


x 2 x 


120 


7 


X 


2 


X 


120 


45149* 


150.15022 


2.475196 


0.692 


spec-z 


2 


x 2 x 60 


2 


x 2 x 


120 


2 


X 


2 


X 


120 


46960 


149.72560 


1.810861 


0.747 


spec-z 


10 


x 2 x 60 


7 


x 2 x 


120 


5 


X 


2 


X 


120 


48730 


150.18233 


1.700827 


0.735 


spec-z 


2 


x 2 x 60 


2 


x 2 x 


120 


2 


X 


2 


X 


120 


48760* 


150.68904 


1.757104 


0.586 


phot-z 


2 


x 2 x 60 


2 


x 2 x 


120 


2 


X 


2 


X 


120 


49150 


150.09402 


2.299131 


0.691 


spec-z 


11 


x 2 x 60 


11 


x 2 x 


120 


5 


X 


2 


X 


120 


49207 


149.96306 


2.431469 


0.663 


phot-z 


5 


x 2 x 60 


8 


x 2 x 


120 


5 


X 


2 


X 


120 


49223 


150.32326 


2.449379 


0.744 


spec-z 


6 


x 2 x 60 


8 


x 2 x 


120 


5 


X 


2 


X 


120 


49946* 


150.63387 


2.593702 


0.658 


spec-z 


2 


x 2 x 60 


2 


x 2 x 


120 


2 


X 


2 


X 


120 


50513* 


149.65569 


2.600808 


0.735 


spec-z 


2 


x 2 x 60 


4 


x 2 x 


120 


3 


X 


2 


X 


120 


50519 


150.35817 


2.638588 


0.703 


spec-z 


10 


x 2 x 60 


11 


x 2 x 


120 


7 


X 


2 


X 


120 


50717 


149.68130 


2.681122 


0.659 


spec-z 


9 


x 2 x 60 


8 


x 2 x 


120 


5 


X 


2 


X 


120 


51182* 


150.62514 


1.802925 


0.626 


spec-z 


2 


x 2 x 60 


2 


x 2 x 


120 


2 


X 


2 


X 


120 


51359 


149.79614 


1.774491 


0.740 


phot-z 


5 


x 2 x 60 


5 


x 2 x 


120 


4 


X 


2 


X 


120 


51415 


149.69719 


1.905216 


0.663 


spec-z 


6 


x 2 x 60 


10 


x 2 x 


120 


6 


X 


2 


X 


120 


51416 


150.20341 


1.902644 


0.753 


spec-z 


4 


x 2 x 60 


2 


x 2 x 


120 


2 


X 


2 


X 


120 


51427 


150.02583 


1.926421 


0.661 


spec-z 


4 


x 2 x 60 


2 


x 2 x 


120 


2 


X 


2 


X 


120 


51429 


149.62653 


1.931427 


0.664 


spec-z 


4 


x 2 x 60 


5 


x 2 x 


120 


4 


X 


2 


X 


120 


51477 


150.48152 


2.011073 


0.660 


spec-z 


6 


x 2 x 60 


7 


x 2 x 


120 


5 


X 


2 


X 


120 


51489 


149.87846 


2.031983 


0.677 


spec-z 


5 


x 2 x 60 


7 


x 2 x 


120 


5 


X 


2 


X 


120 


51497* 


150.44370 


2.049112 


0.668 


spec-z 


4 


x 2 x 60 


7 


x 2 x 


120 


6 


X 


2 


X 


120 


51523 


150.43016 


2.086895 


0.659 


spec-z 


4 


x 2 x 60 


2 


x 2 x 


120 


2 


X 


2 


X 


120 


51538 


149.76747 


2.117411 


0.670 


spec-z 


11 


x 2 x 60 


10 


x 2 x 


120 


5 


X 


2 


X 


120 


51604* 


149.91245 


2.200366 


0.688 


spec-z 


4 


x 2 x 60 


2 


x 2 x 


120 


2 


X 


2 


X 


120 


51645* 


149.57237 


2.262714 


0.764 


phot-z 


2 


x 2 x 60 


2 


x 2 x 


120 


2 


X 


2 


X 


120 


51653 


149.80357 


2.288883 


0.762 


spec-z 


8 


x 2 x 60 


8 


x 2 x 


120 


4 


X 


2 


X 


120 


51728 


150.28973 


2.400017 


0.614 


spec-z 


14 


x 2 x 60 


15 


x 2 x 


120 


9 


X 


2 


X 


120 


51739* 


150.47203 


2.410231 


0.668 


spec-z 


4 


x 2 x 60 


5 


x 2 x 


120 


4 


X 


2 


X 


120 


51852* 


150.05379 


2.589671 


0.697 


spec-z 


2 


x 2 x 60 


2 


x 2 x 


120 


2 


X 


2 


X 


120 


51916 


150.48936 


2.688268 


0.657 


spec-z 


13 


x 2 x 60 


15 


x 2 x 


120 


8 


X 


2 


X 


120 


51982 


150.40459 


2.780556 


0.595 


phot-z 


4 


x 2 x 60 


5 


x 2 x 


120 


4 


X 


2 


X 


120 



0.74 
0.86 
0.89 
0.78 
0.79 
0.79 
1.17 
0.75 
0.85 
0.87 
0.87 
0.81 
0.75 
0.82 
0.81 
0.77 
0.73 
5.02 
1.05 
1.42 
0.86 
0.76 
3.82 
1.23 
1.79 
2.13 
0.87 
1.09 
0.94 
2.16 
1.18 
0.83 
1.19 
3.04 
1.17 
0.84 
2.59 
2.38 
1.22 
1.15 
1.08 
0.88 
2.06 
0.97 
1.84 
7.12 
1.16 
0.71 
1.12 
3.39 
0.76 
1.33 



10.65 
10.54 
10.57 
10.58 
10.72 
10.69 
10.75 
10.52 
10.73 
10.75 
11.00 
10.83 
10.60 
10.53 
10.70 
10.79 
10.60 
11.70 
10.78 
11.03 
10.83 
10.72 
11.40 
11.05 
11.11 
10.94 
10.76 
10.80 
10.95 
11.16 
10.99 
10.72 
10.84 
11.23 
10.98 
10.66 
11.33 
11.19 
10.84 
10.71 
10.67 
10.75 
11.02 
10.60 
11.00 
11.83 
10.78 
10.51 
10.71 
11.38 
10.66 
10.79 



11.06 
11.06 
11.00 
10.97 
10.74 
11.01 
11.15 
10.90 
10.99 
11.06 
10.91 
10.98 
10.97 
11.01 
10.96 
11.12 
10.54 
11.64 
10.39 
11.22 
11.13 
11.04 
11.40 
11.27 
11.29 
10.90 
10.79 
10.73 
10.95 
11.17 
11.00 
11.01 
10.88 
11.35 
10.81 
10.82 
11.32 
11.18 
10.86 
11.18 
11.00 
10.74 
11.24 
10.90 
10.90 
11.59 
11.05 
10.68 
10.93 
11.41 
10.93 
10.82 



0.00 
0.00 
0.00 
0.00 
0.64 
0.00 
0.00 
0.00 
0.00 
0.11 
0.98 
0.00 
0.00 
0.00 
0.00 
0.00 
0.97 
1.00 
0.74 
0.27 
0.00 
0.00 
0.54 
0.14 
0.41 
0.99 
0.25 
0.99 
0.67 
0.76 
0.99 
0.00 
0.17 
0.85 
0.99 
0.09 
0.76 
0.83 
0.52 
0.00 
0.00 
0.99 
0.21 
0.00 
0.53 
1.00 
0.00 
0.00 
0.48 
0.99 
0.00 
1.00 



NOTE. — Column (1): ID in the S-COSMOS MIPS 24 (ira catalog. X-ray sources are indicated by stars. Column (2): Right Ascension 
of the optical counterpart. Column (3): Declination of the optical counterpart. Column (4): Redshift. Column (5): Redshift flag. Columns 
(6)— (8): Exposure times for the 3 IRS modules; the total integration time equals the multiplication of the number of cycles, the number 
of dithering positions of each cycle, and the exposure time of an individual frame (in seconds). Column (9): MIPS 24 fim flux in mJy, 
these are before the color correction described in § 2.1. Columns (10)-(11): Rest-frame luminosities in Log(L0) at 15 fim and IRAC 8 
fim; directly measured from the IRS spectra. Column (12): Contribution of the AGN component to the integrated rest-frame 5.5 to 20 /im 
luminosity. Sources outside of our redshift bin (0.6 < z < 0.8) are included for completeness. 
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TABLE 2 

Star-forming 15/8 /m Luminosity Functions 





W(15^m) 


<f>(l5/j,m) 










8ij)(8(j,m) 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


O) 


(8) 




MIPS Sample 


9.68 


1017 


-2.39 


-3.89 


9.94 


873 


-2.33 


-3.78 


9.82 


850 


-2.50 


-3.97 


10.10 


938 


-2.44 


-3.93 


9.98 


675 


-2.60 


-4.02 


10.24 


757 


-2.55 


-3.99 


10.12 


468 


-2.76 


-4.10 


10.40 


548 


-2.69 


-4.06 


10.27 


283 


-2.98 


-4.21 


10.55 


386 


-2.85 


-4.14 


10.43 


173 


-3.20 


-4.31 


10.69 


224 


-3.08 


-4.26 


10.57 


69 


-3.59 


-4.51 


10.85 


117 


-3.37 


-4.40 


IRS Sample 


10.57 


9 


-3.61 


-4.01 


10.85 


4 


-3.33 


-3.39 


10.73 


14 


-4.10 


-4.67 


10.99 


16 


-3.85 


-4.44 


10.88 


4 


-4.66 


-4.96 


11.15 


6 


-4.47 


-4.85 


11.02 


1 


-5.26 


-5.26 


11.30 


1 


-5.26 


-5.26 



NOTE. — Columns (1,5): Luminosity bin center. Columns (2,6): Number of objects in 
the bin. Columns (3,7): Comoving volume density in log(Mpc~ 3 dex~'). Columns (4,8): 
Poisson error of the comoving volume density in log(Mpc~ 3 dex~' ). 



TABLE 3 

AGN 15/8/iM Luminosity Functions 



L 15 




<f>(l5/j,m) 


<50(15/^m) 


Ls 




0(8/^m) 


84>(Sfj,m) 




10.65 


3 


-4.18 


-4.32 


10.77 


4 


-4.58 


-4.88 


10.80 


6 


-4.44 


-4.83 


10.92 


3 


-4.78 


-5.02 


10.95 


3 


-4.78 


-5.02 


11.07 


3 


-4.78 


-5.02 


11.10 


3 


-4.78 


-5.02 


11.22 


2 


-4.96 


-5.11 


11.25 


2 


-4.96 


-5.11 


11.37 


1 


-5.26 


-5.26 


11.40 


1 


-5.26 


-5.26 


11.52 


1 


-5.26 


-5.26 


11.70 


1 


-5.26 


-5.26 


11.67 


1 


-5.26 


-5.26 



TABLE 4 

Best-Fit Parameters of Star-forming Luminosity Functions 



Parameter 15 fj,m IRAC8/^m 8-1000 fim Unit 



L 


10.31 ±0.04 


10.53 ±0.04 


11.58 ±0.03 


log(L ) 


4>* 


-2.38 ±0.05 


-2.42 ± 0.06 


-2.72 ±0.04 


log(Mpc~ 3 dex" 1 ) 


a 


-1.25 ±0.07 


-1.32±0.07 


-1.46 ±0.05 




a 






0.30 ±0.01 




P [ 


7.66 ±0.02 


7.88 ±0.02 


8.85 ±0.02 


log(L Mpc- 3 ) 



NOTE. — Total integrated luminosity from the analytic fit. 



TABLE 5 
8-1 000 fiM Luminosity 
Function 





N 





50 




11.60 


10 


-2.83 


-3.21 


11.75 


44 


-3.08 


-3.86 


11.90 


46 


-3.44 


-4.23 


12.05 


27 


-3.78 


-4.43 


12.20 


11 


-4.38 


-4.90 


12.35 


6 


-4.66 


-5.04 


12.50 


1 


-5.43 


-5.43 


12.80 


1 


-5.43 


-5.43 



